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Abstract. Wo evaluate the chemical potential of a one-dimensional quantum 
dot, coupled to an atomic force microscope tip. The dot is described within the 
Luttingcr liquid framework and the conductance peaks positions as a function of 
the tip location are calculated in the linear and non-linear transport regimes 
for an arbitrary number of particles. The differences between the chemical 
potential oscillations induced by Friedel and Wigner terms are carefully analyzed 
in the whole range of interaction strength. It is shown that Friedel oscillations, 
differently from the Wigner ones, are sensitive probes to detect excited spin states 
and collective spin density waves involved in the transport. 
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When electrons are confined in a finite portion of condensed matter, such as in 
quantum dots [JJ, many intriguing physical effects appear. Quantum confinement 
induces Friedel oscillations [2 in the electron density with a typical wavelength 
(2kp)~ 1 (fcp the Fermi momentum). In addition, when the Coulomb repulsion 
dominates over the kinetic energy, electrons tend to form a Wigner crystal [31 0] giving 
rise to density oscillations at wavelength (4&f) . In two-dimensional dots, correlation 
effects have been extensively studied mainly by means of numerical techniques [3J HJ 
HJini 13 • Due to the reduced dimensionality, one-dimensional (ID) quantum dots, such 
as carbon nanotubes [S], show even more dramatic interaction effects. They have been 
indeed the subject of intense numerical investigations [H EH EH E2 EH ED- Despite 
their precision, numerical methods are usually restricted to a low number of particles 
due to their complexity. An analytical approach, widely employed to describe the 
physics of interacting ID electrons, is the Luttinger liquid model [T51 UM ITTl IT81 [T5] . 
In connection with the bosonization technique, it represents a powerful method to 
study the low-energy physics, allowing to explore the limit of large particle numbers 
and addressing their transport properties [201 EU HH EH1 EI] and the formation of 
Wigner crystals ESI ESI E3 EH E2] . 

One of the key tools to investigate the emergence of Friedel and Wigner 
correlations is exploiting the dot transport properties in the presence of a local 
probe. Scanning tunnel microscope tips have been proposed in the past to detect 
Friedel [3U1 EJJ GUI or Wigner oscillations [33] and even local electron-vibron 
coupling [34 . The most natural and powerful choice is however that of an atomic force 
microscope (AFM) tip which capacitively couples to the dot density and thus allows to 
capture its oscillations, a technique called "Scanning gate microscopy" [351 ISS]- Recent 
theoretical studies ESI [37] have investigated the influence of an AFM tip on the dot 
chemical potential revealing its spatial oscillations as a function of the tip position. 
These works, however, have only considered the linear regime, allowing to access only 
the ground state properties of systems with low electrons numbers [13] . 

In this paper we will go beyond these limitations. Employing the Luttinger liquid 
and the bosonization language ESI EH ESI EHl I2S] we will develop a fully analytical 
framework which will enable us to study the position of the conductance peaks in the 
presence of a moving AFM tip, both in the linear and in the non-linear regime. To this 
end, we will evaluate the chemical potential traces as a function of the tip position and 
electron interaction strength, without any limitation on the electron number. We will 
show that Friedel and Wigner oscillations produce different signatures due to the spin 
dynamics, demonstrating that Friedel correlations are sensitive to the spin populations 
of electrons in the quantum dot. In the non-linear case, where higher excited spin are 
involved, Friedel oscillations increase their number, reaching the one of the Wigner 
case for a fully polarized dot. Also collective spin density excitations induce additional 
density modulations which can as well be detected in the chemical potential traces. 

We consider an interacting ID quantum dot, treated as a Luttinger liquid with 
open boundaries. The reference state is set with an even electron number N° and 
Fermi momentum kp = ttN°/2L (L dot length). The low energy Hamiltonian takes 
the form [18 (H = 1) 



Here, d v>q are the boson operators of the collective charge (v — p) and spin [y = a) 
density waves, with quantized wave number q = im q /L, n q £ N* and energy e v {q) = 




(1) 



f=p,<y q>0 
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e v n q , where e v — ttv v /L. They propagate at a renormalized velocity v v — Vp/g v , (vf 
the Fermi velocity) without a change in the total number of charges and spin. The 
repulsive electron interaction strength is represented by g p = g < 1, while g a = 1 
for an SU(2) invariant theory [IB]. The second part of the Hamiltonian represents 
the zero mode contributions with AN PtCr = AN + ± A7V_, where AN S = N s — N® 
are the extra electrons per spin direction (s = ±) with respect to the reference 
Nl = N°_ = N°/2. The total charge and spin are N p = AN p + N and N a = AN a , 
with energies E p and E a . Despite the microscopic model provides their quantitative 
estimates {E u = TTv v j2Lg v ), several effects, e.g. coupling with the external gates or 
long range interactions, cause deviations especially in the charge sector. Therefore, 
we treat E p as a free parameter, with E p ^> E a , while keeping E„ = ttvf/2L = e a /2. 
Let us denote the eigenstates of the above Hamiltonian as \S) = \N p , N a , {n?}, {n q }} 
with energies 



m = E 



-AN?, 



E 

q>0 



(2) 



where {n q } are the occupation numbers for the collective mode v at momentum q. 
The electron operator fy s (x), with boundary conditions 'I's(O) = ^ S (L) = 0, 
is represented in terms of right-moving electrons as ^> s (x) — e lkFX ip Si R(x) — 



In bosonized form, 



e L e 



Here, {rj s ,r) s /} 



a = hp 1 the cutoff length, [9 S , AN s i] = iS s s > and 



<m*)=E 



q> 



,-aq/2 



[(cos (qx) — ig v sin (qx)) + h.c] 



(3) 



(4) 



The total electron density p(x) = J2 S Ps( x ) with p s (x) = ^l(x)^ s (x) consists of 
several contributions [35] HI] ■ Here we take into account the most relevant adopting a 
number-conserving formalism |29j for open boundaries. The smooth long-wave term 



H s Ps w ( x ) is g iven b y 



iff 



(x) 



hp 
n 



+ -d x tp s (x) - —d x h(x) 

L TT TT 



with ip s (x) = \f p (x) + s<p a (x)] /V2 where 



(5) 
(6) 



sin(27ra; / L) 



and h(x) = I tan 1 ; ; — 

The oscillating Friedel contribution is p F (x) oc J2 s Ps( x ) with Px( x ) 



rs + h.c. In bosonic representation one has 



where 



pf( x ) 



2tt 



d x sm[£(AN s ,x,g) + 2ip s (x)} , 



£(n, x, g) = 2kpx H 2g h(x) . 



(7) 



(8) 



In addition to these "standard" terms we include the so-called Wigner contribution 

p w (x) cx e~ AlkpX ^ + R (x)ip + ^R(—x)ipl_ r {x)iP-^r{— x) + h.c. which may arise due to 
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interaction effects, band curvature or other external perturbations [26, 28, 27 . In the 
bosonization language [25] one has 



p w (x) 



1 

"2tt 



d x sin 



2C 



AN n 



(9) 



Note that these terms constitute the most relevant contributions to the electron 
density, whose amplitudes are to be interpreted as model parameters (23 |2H] ■ The 
total density, satisfying boundary conditions p(0) — p(L) — 0, can then be expressed 
as 

p(x) = p LW {x) + (1 - \)p F (x) + \p w (x) , (10) 

with A €E [0, 1] a free parameter. 

We now analyze the expectation values p GS (x) — (S\p^(x)\S)cs on the ground 
state \S) G s = \N^ s ,N^ s ,{ n P = 0},{n° = 0}) for the different contributions 
£ = LW, F, W. Note that, due to N° being even, one has N^ s = N° and N^ s = 
with AN^ S = for an even number of electrons, while N^ s = N° + 1 and N^ s = ±1 
with AN^ S — 1 for the nearest larger odd electron number. Using the bosonization 
technique one obtains 

2k F AN° S 



Pgs( x ) 



- 2 y -d x h{x) . 

L 7T 



,9)} 



Pgs( x ) = 



- —djK w sin 

Z7T 



2C 



(ANf s 



,x,g 



with the enveloping functions (rjp = 



_ i+fl . 

2 ' 



sinh 



7ra 
2L , 



Ww = 2g) 



^sinh 2 (ff)+sin 2 (^) 



(12) 
(13) 

(14) 



Let us now study the Friedel and Wigner contributions separately. Figure [T] shows 
the electron density, see Eq. ( 10 ), in the ground state with the long wave part and the 
Friedel (A = 0) or the Wigner (A = 1) terms. In both cases it exhibits an oscillatory 
behaviour but with different patterns and amplitudes. The Friedel contribution shows 
oscillations related to the two different electron spin populations. In particular, for 
even electrons one has N^ s = N^ s /2 and thus the superposition results in N^ s /2 
peaks. For odd numbers, one of the two subpopulations has (N^ s + l)/2 electrons, 
while the other has (N^ s — l)/2. The superposition for both spin directions has then 
(N^ s + l)/2 peaks. The Wi gner correction, on the other hand, is insensitive to spin, 



see Eq. (13), and depends on N^ s only, which is also the number of observed peaks. 



Concerning amplitudes one can easily see that Friedel oscillations are dominant over 
the Wigner ones for weak interactions, Panels (a,c), while for strong interactions, 
Panels (b,d), the situation is reversed. This fact is particularly evident when the 
number of particles increases. 

Let us now consider the combined effects of both contributions. Figure [2] shows Pgs( x ) 
as a function of x and g for A = 0.5, see Eq. (10). For weak interactions (g s» 1) 
one can distinguish (N^ s + l)/2 peaks (in both panels N^ s is odd). At strong 
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Figure 1. Ground state dot electron density (units L ') as a function of x (units 



L) for (a) 



5, g = 0.9; (b) N% s = 5, g = 0.1; (c) N^ s = 20, g 



0.9; 



(d) K 



20, g = 0.1. The blue (red) curve represents the Friedel (Wigner) 



contributions in addition to the long-wave part. In (a,b) a = 2n/L, in (c,d) 
a- 1 = IOtt/L. 




interactions (g — > 0), Wigner correlations grow and eventually the density exhibits 
sharp oscillations with N^ s peaks. Thus, even when Wigner correlations are mixed 
along with Friedel ones, the presence of strong interactions makes the Wigner channel 
the relevant one. The abov e b ehaviour can be attributed to the enveloping functions 
K F (x) and K w (x) in Eq. (14). For g — >• 1, the weaker power law in K F (x) (namely, 
^-(i+s)/ 2 ) m comparison to that of K w (x) (namely, x~ 2g ) leads to a suppression 
of the Wigner channel. On the other hand, when g — > one has K w (x) — > 1 and 
K F (x) <C 1 away from the borders. Thus, for strong interactions Friedel oscillations 
are damped, while the Wigner ones are still fully developed. 

We would like now to show how to probe electron density correlations in the 
presence of a movable AFM tip. This has already been considered as a powerful 



Non-linear Coulomb blockade microscopy of a correlated one- dimensional quantum dot6 



tool to investigate the electronic correlated state for linear transport in dots with few 
electrons |13j . Here, we will consider the more general non-linear regime, without any 
constraint on the number of electrons. 

The charged AFM tip is capacitively coupled to the dot at position x, with coupling 
H t ip = Vop(x), assuming a tip narrower than the average wavelength of the density 
oscillations. In the sequential regime, for a given transition \S) — > \S'} electrons tunnel 
between the dot and the leads and the number of dot charges oscillates between N p 
and N p + 1, with the spin constraint \N' a — jV^I = 1. The onset of a transition is 
signalled by peaks in the differential conductance. In the following we will determine 
their positions as a function of the tip. The key quantity to evaluate is the generalized 
chemical potential 

Vs^s> (x) = E tot (S',x) - E tot (S,x) , (15) 

defined in terms of the total energy E to t(S,x) in the presence of the tip. Indeed, 
whenever (J,s^S' (x) ~ P"xi where ^ s the electrochemical potential of the lead x 
the given transition becomes allowed. We will evaluate Eq. (15) in the weak tip- 
coupling regime, namely E tot (S, x) — E(S) + SE(S, x) with SE(S, x) the lowest order 
correction. 

We will consider only non-degenerate energy levels, whose corrections to the energy 
E(S) in Eq. ([2| are given by SE(S,x) = Vo(S\p(x)\S) . This is a relevant case since 
it applies to zero- mode spin excited states (the degeneracy on ±|iV CT | is explicitly 
diagonal and does not cause problems) and to the lowest lying spin density waves with 
a singly occupied bosonic state. The latter is indeed non degenerate in the presence 
of interactions since e a < s p . Generalizing the results obtained for the ground state 
one has 5E(S, x) = 5E LW (S, x) + (1 - X)SE F (S,x) + \6E w (S,x), with 



5E LW {S,x) = Vi 



SE F {S,x) 




ANp 
L 



2—d x h{x) 

7T 



K F G F J2^[£(AN s ,x,g)] 



6E w (S,x) 



Vq 

'2tt 



s=± 



■d x [K w G w sm 



x, g 



(16) 



(17) 



(18) 



Here we introduced G F (x) = JX ll g >o L K [A^xf] and G w (x) = U q>0 K% [A p q {xf] , 
with L n (x) the Laguerre polynomials stemming from generalized Franck-Condon fac- 
tors [301131] and A^(x) ~ \J^j^e~ aq ^ 2 sin(gx). In the absence of collective modes 
one has G F (x) = G w \x) = 1. 

The corresponding chemical potential is then decomposed as the bare one, fJ-o.s^S' = 
E(S')-E(S) - see Eq. ^ - and the tip correction, 5p, s ^S'(x) = 5E(S',x)-SE(S,x) 
- see Eqs. ( 16fl8 1. In short notation, omitting the states, one has 

li(x) =no + 5p,(x) (19) 

(20) 



5fi{x) = 5fi LW + (1 - A) V (a) + \5n w {x) 



w i 



We are now in the position to investigate the chemical potential variation as 
a function of the tip position for a given transition. We start by considering the 
linear regime where only the ground states are involved: |5)gs with total number N p 



and 2|7V ct | = i _ (-1) w p, and \S') GS with NL = N p + 1 and 2|iV;| = 1 - 
Figure [3] shows the correction 5/j,(x) in Eq. (20). One can observe sharp oscillations, 
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Figure 3. Chemical potential corrections <5/i(ie) (units Vq/ L) in the linear regime 
as a function of x (units L) for the transition with N p = 16. (a) Friedel corrections 
only (A = 0) and g = 0.7 ; (b) same as in (a) but for Wigner (A = 1); (c) same as 
in (a) but with Friedel and Wigner (A = 1/2); (d) same as in (c) but for stronger 
interactions g = 0.1. In all panels, a~ 1 = 8n/L. 



induced by the difference of two oscillating quantities and in contrast to the weaker 
ones displayed in the electron density. The Friedel term, shown in Panel (a) for even 
N p and N a = 0, displays 1 + N p /2 maxima and N p /2 minima. The number of these 
maxima is related to the spin population of the dot final state. Indeed, the initial 
spin populations are N s = N p /2. When an extra electron, with spin up, enters into 
the dot, the state with N' a = 1 is reached with N+ = 1 + N p /2 and N'_ = N p /2. 
These are indeed the values which correspond to the number of maxima and minima 
exhibited. A similar behaviour holds for an odd number N p (not shown) where one 
finds (N p + l)/2 maxima and (N p — l)/2 minima. Thus, in the Friedel case the 
oscillations of the chemical potential are a sensitive probe to detect the spin sub- 
populations of the electron involved in the linear transport process. The Wigner case, 
on the other hand, is not sensitive to the spin direction. As shown in Panel (b), it 
always displays N p + 1 maxima and N p minima, allowing to count the total number of 
electrons only. Panels (c,d) show the combined effect of equally weighted Friedel and 
Wigner corrections. For strong interactions, Panel (d), Wigner corrections are clearly 
well pronounced, differently from the weak interactions case, Panel (c). These results 
are in agreement with the ones obtained for few particles with an exact diagonalization 
procedure [13] . 

Let us turn to the non-linear case, with excited zero mode spin states and for 
the moment without collective modes. This has never been investigated in literature 
to our knowledge. The transitions are between states of the form |<S) = \N p , N^} and 
\S') = \N p + l,N' a ), with N a , N' a not restricted to the ground state values anymore. 
We select, in particular N' a = A^ + l and N p = 2k with n 6 N*, focusing on the Friedel 
contribution only. The Wigner one is indeed insensitive to the total spin - see Eq. ( 18 ) - 
and will give, for each of the above transitions, a contribution equal to that of the linear 
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case. The possible spin transitions are |iV p ,0) — > \N P + 1,1); \N P + 1,1) — > \N p ,2); 
• • • ; \N p , N p ) —> \N p + 1, N„ + 1). Namely, the dot starts in the ground state with 
N s = Np/2, eventually becoming fully polarized. As a representive example we show 
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Figure 4. Chemical potential /i(x) (units E a ) as a function of x (units L) 
including only the Friedel corrections, (a) Traces of the transition |6,0) — > |7, 1} 
(blue), |6,2) -» [7,3) (green), |6,4) -* |7,5) (orange), |6,6> -¥ |7,7) (red); (b) 
Traces for a transition |6, 0,{n p q = 0},{n£ = 0}> -> |7, l,{n p q = 0},{nj = 
Sq.n/L}) involving a spin density wave (black) and |6,2) — > |7, 3} involving zero- 
modes only (green). In all panels g = 0.7, E p = 10-E CT , Vb = 0.08E a and 
a - 1 =3n/L. 



in Fig. |4|a) the position dependence of the chemical potential in Eq. ( 19 1 with the 
Friedel contribution for the tunnel-in transitions. Starting with N p /2 + 1 maxima in 
the ground states (blue trace), an increasing number of oscillations appear as higher 
spin states get involved. Eventually, the transition involving fully polarized dot state 
is reached and the chemical potential shows with N p + 1 maxima (red trace). Note 
that this is the same number of maxima that one would obtain from the Wigner 
corrections. 

The sensitivity of the Friedel corrections to spin suggests to employ non-linear 
transport experiments in the presence of an AFM tip to explore highly excited 
spin states of quantum dots. This seems particularly desirable for not too strong 
interactions, when the Wigner corrections do not dominate over the Friedel ones. It is 
also possible to conclude that the observation of N p + 1 peaks in the chemical potential 
for an excited transition does not provide alone a clear-cut evidence of the presence 
of Wigner correlations. 

We close our discussion briefly addressing how the above effect can be obtained 
also for transitions involving excited spin density waves previously neglected. In 
general, and in contrast to zero-mode excitations, collective modes are assumed to 
relax to the ground state over a time shorter than the average tunneling time of 
electrons due to several possible types of dissipative coupling. This implies excitations 
only in final states. To be more specific we consider a transition from the ground state 
to the lowest-lying collective spin excitation with q = tt/L. The two involved states 
are \N P , 0, {0}„ {0} CT ) -> \N p + 1, 1, {0} p , {n° q = 5 q ^ /L }) with N p even. Figure JJb) 
shows the corresponding chemical potential with Friedel corrections for N p = 67 It 
exhibits five peaks and the overall behaviour is similar to the green trace in Fig. [4^a.) 
corresponding to |6, 2) — > |7, 3). Indeed since 2E a = e a the two transitions have the 
same bare chemical potential (j,q = E p (2N p + 1) /2 + 5E a /2 but different corrections 
5fi(x) also due to the different contribution of the function G F (x) in Eq. (17). 

In conclusion, we have studied the chemical potential traces of an interacting ID 
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quantum dot in the presence of a coupling with a moving AFM tip. As a general 
trend, we have found that Friedel modulations dominate for weak interactions and 
are sensitive to spin populations, while the Wigner oscillations become relevant at 
strong interactions and depend on the total charge sector only. We demonstrated that 
this results in markedly different behaviours. In the linear regime, Friedel correlations 
exhibit half of the number of oscillations than the Wigner ones. On the other hand, in 
the non-linear case, when higher excited spin are involved, they increase the number of 
oscillations, reaching the one of the Wigner case for a fully polarized dot. We expect 
that predicted oscillations of the chemical potential traces as a function of the AFM 
tip can be observed in transport experiments. 
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